Accuracy of the post-Newtonian approximation. II. Optimal asymptotic expansion of 
the energy flux for quasicircular, extreme mass-ratio inspirals into a Kerr black hole. 
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We study the effect of black hole spin on the accuracy of the post-Newtonian approximation. We 
focus on the gravitational energy flux for the quasicircular, equatorial, extreme mass-ratio inspiral 
of a compact object into a Kerr black hole of mass M and spin J. For a given dimensionless spin 
a = J /M^ (in geometrical units G = c = 1), the energy flux depends only on the orbital velocity 
V or (equivalently) on the Boyer-Lindquist orbital radius r. We investigate the formal region of 
validity of the Taylor post-Newtonian expansion of the energy flux (which is known up to order 
beyond the quadrupole formula), generalizing previous work by two of us. The error function 
used to determine the region of validity of the post-Newtonian expansion can have two qualitatively 
different kinds of behavior, and we deal with these two cases separately. We find that, at any fixed 
post-Newtonian order, the edge of the region of validity (as measured by u/fisco, where Hisco is 
the orbital velocity at the innermost stable circular orbit) is only weakly dependent on a. Unlike 
in the nonspinning case, the lack of sufficiently high order terms does not allow us to determine 
if there is a convergent to divergent transition at order . Independently of a, the inclusion of 
angular multipoles up to and including ^ = 5 in the numerical flux is necessary to achieve the level 
of accuracy of the best-known {N = 8) PN expansion of the energy flux. 



I. INTRODUCTION 

Binaries of compact objects, such as black holes (BHs) 
and/or neutron stars, are one of the main targets for 
gravitational-wave (GW) observations. When the binary 
members are widely separated, their slow inspiral can be 
well-described by the post-Newtonian (PN) approxima- 
tion, a perturbative asymptotic expansion of the "true" 
solution of the Einstein equations. The small expansion 
parameter in the PN approximation is u/c, where v is 
the orbital velocity of the binary and c is the speed of 
light. Asymptotic expansions, however, must be used 
with care, as the inclusion of higher-order terms does not 
necessarily lead to an increase in accuracy. Therefore one 
would like to determine the optimal order of expansion 
and the formal region of validity of the PN asymptotic 
series [H i-e. the order and region inside which the 
addition of higher order terms increases the accuracy of 
the approximation in a convergent fashion. 

In a previous paper (Q, henceforth Paper I), two of 
us investigated the accuracy of the PN approximation 
for quasicircular, nonspinning (Schwarzschild), extreme 
mass-ratio inspirals (EMRIs). By comparing the PN ex- 
pansion of the energy flux to numerical calculations in the 
perturbative Teukolsky formalism, we concluded that (i) 
the region of validity of the PN expansion is largest at 
relative 3PN order - i.e., order {v/cf' (throughout this 
paper, a term of 0{v'^^) is said to be of A^PN order); 
and (ii) the inclusion of higher multipoles in numerical 
calculations is necessary to improve the agreement with 
PN expansions at large orbital velocities. The fact that 
the region of validity is largest at 3PN could be a hint 
that the series actually diverges beyond 3PN order, at 



least in the extreme mass-ratio limit. 

This paper extends our study to EMRIs for which the 
more massive component is a rotating (Kerr) BH. The 
present analysis focuses on the effect of the BH spin on 
the accuracy of the PN expansion. We generalize the 
methods presented in Paper I to take into account cer- 
tain pathological behaviors of the error function, used to 
determine the region of validity. This generalization may 
also be applicable to comparable-mass systems. 

A surprising result we find is that the edge of the region 
of validity (the maximum velocity beyond which higher- 
order terms in the series cannot be neglected), normalized 
to the velocity at the innermost stable circular orbit, is 
weakly dependent on the Kerr spin parameter a. In fact, 
this edge is roughly in the range w/wisco S [0.3,0.6] for 
almost all PN orders, irrespective of a. This suggests, 
perhaps, that the ratio ?;/t;isco is a better PN expansion 
parameter than u/c, when considering spinning BHs. 

Another surprising result is related to the behavior of 
the edge of the region of validity as a function of PN 
order. In the nonspinning case, two of us found that 
beyond 3PN order, 0{v^ /c^), this edge seemed to consis- 
tently decrease with PN order [Ij. This was studied up 
to ©(w^^/c^^), the largest PN order known for the non- 
spinning case. In the spinning case, however, the series 
is known only up to 0(u^/c^), and we are thus unable 
to conclusively determine if the trend found in the non- 
spinning case persists. Higher-order calculations will be 
necessary to draw more definite conclusions. 

Numerical (or in this case, perturbative) calculations 
of the energy flux rely on multipolar decompositions of 
the angular dependence of the radiation. By comparing 
the convergence of the multipolar decomposition to the 
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convergence of the PN expansion of the energy flux, we 
find that for w/c ~ 0.1 the inclusion of multipoles up to 
and including i = 5 seems necessary to achieve the level 
of accuracy of the best-known {N = 8) PN expansion of 
the flux. These conclusions are also largely independent 
of the spin parameter a. 

The rest of this paper is organized as follows. In Sec- 
tion |lT] we present the energy flux radiated by quasicircu- 
lar, equatorial Kerr EMRIs in the adiabatic approxima- 
tion, as computed in PN theory and with accurate 
frequency-domain codes in BH perturbation theory ■ 
In Section Uni we discuss the region of validity of the PN 
approximation in terms of the normalized orbital veloc- 
ity t;/uisco and of the normalized orbital radius r/risco, 
where ISCO stands for the innermost stable circular or- 
bit. We consider both corotating and counterrotating 
orbits. In Section ITVl we study the number of multipolar 
components that must be included in the numerical flux 
in order to achieve sufficient accuracy. Finally, in Section 
|V]we present our conclusions. We follow the same nota- 
tion as in Paper I. In particular, from now on we will use 
geometrical units {G = c— 1). 



II. ENERGY FLUX FOR QUASICIRCULAR, 
EQUATORIAL EMRIS IN KERR: NUMERICAL 
AND PN RESULTS 

In the PN approximation, the GW energy flux radiated 
to infinity by a test particle in a circular orbit and on the 
equatorial plane of a Kerr BH is given by 0411 
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(1) 



This flux is known up to = 8 when including spins, 
and up to iV = 11 in the nonspinning case. The leading 
(Newtonian) contributiorQ is 



(2) 



where fj, and M are the test particle mass and Kerr BH 
mass, respectively. As we are here interested in the ac- 
curacy of the PN approximation, we will ignore the flux 
of energy going into the horizon, which cannot always be 
neglected when building waveform templates. 

The expansion coefficients Ofc and bk contain both spin- 
independent and spin-dependent terms, where the dimen- 
sionless spin parameter a is related to the Kerr BH spin 
angular momentum via J = aM"^ . These coefficients can 
be found in Eq. (G19) of H, so we do not list them explic- 
itljQ. Note that logarithmic terms only appear at 3PN 



^ Notice that there is a typo in Eq. (18) of Paper I. 
^ See also their Eq. (3.40), that provides a similar expansion in 
terms of the PN orbital velocity parameter v' = \/M/ro. 
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FIG. 1. ISCO velocity (disco, top panel) and radius 
{risco/M, bottom panel) as a function of a. Here and else- 
where we use the convention that a negative spin parameter 
refers to counterrotating orbits. 



and 4PN (i.e., be ^ and bs ^ 0), and that the {ak,bk) 
for 8 < fc < 11 are known only in the spin- independent 
limit. 

Throughout this paper v is the orbital velocity, defined 
as u = (Mfi)^/"^ (where Q is the small body's orbital 
frequency), and related to the Boyer-Lindquist radius r 

by 



r_ _ (1 - at;3)2/3 
M ~ 



whose inverse is 



(r/A/) 

At the ISCO we have [ll| 

''isc 



3/2 



-1/3 



M 



= Z + Z2- v/(3-^i)(3 + Zi + 2Z2), 



Zi = l + (l-a2)i/3 (i + a)i/3 + (l-a)i/3 
Z2^(3a2 + z2)i/2, 



(3) 



(4) 



(5) 
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where a > (a < 0) corresponds to corotating (counter- 
rotating) orbits. 

Using Eqs. ([3]) and ([5]), we also have 
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,,3N2/3 



3 + Z2 T V(3-^i)(3 + ^1+2^2) 



The ISCO velocity can be found by replacing 

^isco in 

Eq. ^ into the velocity-radius relation ([4]). The velocity 
Visco and the radius risco/M are displayed graphically 
in Fig. [1] Observe that, although r^sco M as a ^ 1, 
Visco is bounded by 2'^/^ ~ 0.79. 

The rigorous definition of velocity is a tricky business 
in general relativity. We have here chosen to define veloc- 
ity in a quasi-Newtonian fashion, in terms of the angular 
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FIG. 2. Gravitational energy flux (normalized to Fncwi) as a function of the normalized orbital velocity, ii/iJisco- The left 
panel is for corotating orbits, and the right panel for counterrotating orbits. Different insets refer to different spin parameters 
a, as indicated. The thick black line is the numerical flux. Other linestyles refer to different PN approximations: (thin 
black), F(^) (long-dashed red), F^*'^ (dash-dotted green), F'^^ (dash-dash-dotted blue), F^^^ (dash-dot-dotted orange), F^^' 
(dotted dark green), F'*^ (short-dashed violet). 



velocity and Kepler's \a.w. One can think of this velocity 
as that which would be measured by an observer at spa- 
tial infinity. On the other hand, one can also study the 
velocity measured by an observer in the neighborhood 
of the BH and that is rotating with the geometry; this 
quantity would differ from v in Eq. Q, and, in fact, its 
associated Wisco would tend to 1/2 in the limit a — > 1 (see 
eg. Eq. (3.116) in [ni). This shows that the a ^ 1 limit 
is very delicate, and the precise value of the velocity field 
is an observer-dependent (and a non-invariant) quantity. 
However, once a definition is chosen, the velocity is a per- 
fectly good quantity to parametrize the structure of the 
PN series. 

A first guess at the asymptotic behavior of this series 
can be obtained by simply plotting different PN approx- 
imants F^^^^ and comparing them with high-accuracy, 
numerical results for the energy flux, obtained from a 
frequency-domain Teukolsky code (see [1, for early 
work in the Schwarzschild case, and Fig. 9 in [l^ for 
a related discussion in the Kerr context). The numeri- 
cal results used in this comparison are the same as those 
used in Refs. to study the accuracy of a resummed 
effective-one-body version of the PN approximation to 
model EMRIs. They consist of numerical fluxes, evalu- 
ated for spin parameters ranging from a = 0toa = 0.9 
in steps of Aa = 0.1 (in fact, we also have access to the 
counterrotating flux for a = —0.99). The typical accu- 
racy of these fluxes is better than one part in 10^° for all 
velocities and spins. We refer the reader to Section IIB 
of Q for a more detailed description of the code. 

Figure [2] compares the different PN approximations to 
the numerical flux. The left panel refers to corotating or- 
bits, and the right panel to counterrotating orbits. Differ- 
ent insets correspond to different values of the BH spin, 
and different linestyles represent different orders in the 



PN expansion. As stated earlier, in this figure and in 
the rest of this paper, we neglect energy absorption by 
the BH. Observe that, as first noted by Poisson in the 
Schwarzschild case [1] , the behavior of the PN expansion 
is quite erratic. For any given a, rather than converging 
monotonically, higher-order approximations keep under- 
shooting and overshooting with respect to the "exact" 
numerical result. This oscillatory behavior is quite typ- 
ical of asymptotic expansions, and it has been studied 
in depth, especially for extreme mass-ratio inspirals into 
nonrotating BHs [3, Q. Various authors proposed dif- 
ferent schemes to accelerate the convergence of the PN 
expansion, including Padc resummations p^ - [T^ and the 
use of Chebyshev polynomials [l^ . The asymptotic prop- 
erties of these resummation techniques are an interesting 
topic for future study. 

This figure provides some clues about the edge of the 
region of validity of the PN approximation. For coro- 
tating orbits (left panel of Fig. [5]), as the spin increases 
from zero to a = 0.9, the higher-order PN approximants 
start to deviate from numerical results at lower values 
of w/wisco: this happens roughly when w/wisco — 0-6 for 
a = 0, and when w/wisco — 0.4 for a = 0.9. This leads us 
to naively expect a shrinking of the region of validity of 
the PN approximation as a function of positive a. This 
expectation will be validated (at least qualitatively) in 
Section Hm cf. the bottom-right panel of Fig. U] below. 

At first sight, the results for counterrotating orbits 
(right panel of Fig. [5]) seem surprisingly good. In par- 
ticular, the 3PN approximation (dash-dot-dotted, orange 
line) is almost indistinguishable from the numerical re- 
sult all the way up to w = Wisco when the spin is large. 
Such a good performance is simply because of the well- 
known, monotonically-increasing behavior of iiisco with 
spin, with a minimum as a ^ — 1 (cf. Fig. [T|). Since 
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FIG. 3. Left: absolute value of the remainder of the iV = 3 PN flux, \F ~ F^''^^/ Fucv^t (solid line), and the iV = 4 term 
_ F''^' |/_FNewt (dashed red line). The inset shows the modulus of their difference, Eq. ((H]), in a semilogarithmic scale. 
Right: same as the left panel, but for the N = 6 remainder and N = 7 term. All curves in this plot refer to the counterrotating 
case with spin a = —0.5. The lower (vi, more conservative) and upper {v, less conservative) edges of the region of validity are 
(somewhat conventionally) delimited by the vertical lines, as explained in the main text. 



counterrotating orbits probe a smaller range in v/c (up 
to v/c ^ 0.35 for fast-spinning BHs), the PN approxi- 
mation is more accurate. Unfortunately, prograde accre- 
tion is likely to be more common than retrograde accre- 
tion in astrophysical settings (see e.g. (3). Moreover, 
the 3.5PN and 4PN approximants are significantly worse 
than the 3PN one at a = —0.99. This is consistent with 
the PN series being an asymptotic expansion, as one of 
the characteristic features of the latter is that beyond a 
certain optimal order, higher-order approximants become 
less accurate 



III. REGION OF VALIDITY 

Let us now turn to determining the region of validity 
of the PN approximation for different values of the BH 
spin. For a complete review of asyrnptotic approximation 
techniques wc refer the reader to Paper I presents a 
short introduction to the topic in the present context. 
As explained in those references, the edge of the region 
of validity is determined by the approximate condition 



(7) 



where F denotes the "true" (numerical) result for the 
GW energy flux and F*^^' denotes the A^-th order PN 
approximation. 

An inherent and intrinsic ambiguity is contained in 
Eq. ([7]), encoded in the order symbol. This makes any 
definition of the region of validity of an asymptotic series 
somewhat imprecise. As shown in Fig. [3] (or in Fig. 8 of 
Paper I), there are two qualitatively different scenarios: 

i) Left panel of Fig. O The next-order term — 
starts off smaller than the remainder \F — 



F^^^l, but eventually they cross and separate. We 
can then estimate the edge of the region of validity 
V by solving S^^\v) — 0, where 



(8) 



is the error function. If we also define a more con- 
servative lower edge of the region of validity, vi , as 
the point where 



d,5(^)(i;) 



dv 



(9) 



we can then introduce an uncertainty width of the 
region of validity: Sv = v — vi; see the inset of the 
left panel of Fig. |31 

ii) Right panel of Fig.[3l The remainder and the next- 
order term are of the same order for sufRciently low 
velocities, until eventually the curves separate for 
larger velocities. This situation is the rule, rather 
than the exception, for the problem we consider in 
this paper. When this happens, method i) cannot 
be applied, because (5^^^(w) = has no solutions. 
Given the approximate nature of the order relation- 
ship in Eq. ([7]), we can define the region of validity 
as the point v such that d^'^^v) = Sq, where 5q is 
some given tolerance defined below. 

Higher-order approximations should be sensitive to a 
smaller tolerance, which implies that Sq cannot be set 
arbitrarily. Instead, Sq should be given by the error in the 
difference between the A^th remainder and the (A^-l-l)th- 
order term. This error is presumably of the order of the 
error in the (A'^-l-l)th-order term, and it can be estimated 
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FIG. 4. Edge of the region of validity as a function of a for different PN orders, in the corotating (black, straight line) and 
counterrotating (red, dashed line) cases. The blue, dashed lines for A'^ = 3 refer to the counterrotating case, and they were 
obtained by an alternative method (see the discussion around Fig. [7| below). 



by the {N + 2)th-ordcr term. The imprecision of the 
order symbol is now encoded in the fact that 5q depends 
on V. We can try to estimate its value by evaluating the 
(iV -I- 2)th order term in the middle of the allowed range, 
that is, at t;isco/2: 



So 



F("+2)(^;:sco/2)-F(^+l)(«Isco/2) 



(10) 



This estimate of 5q is not exact, so we can try to provide 
a more conservative lower edge of the region of validity, 
vi, by imposing the conditiorQ 5'^^\vi) = do/2. We can 
then define an uncertainty on the region of validity Sv = 
\v — vi\. This is illustrated pictorially by the vertical lines 
in the right panel of Fig. [31 

Let us now discuss the behavior of the edge of the 
region of validity as a function of the PN order N and 
of the BH spin a. The corotating and counterrotating 
regions of validity and the associated errors are shown 
in Fig. U with solid black (dashed red) error bars for 
corotating (counterrotating) orbits, respectively. 

Let us first consider the corotating case (solid black 
error bars). All results were obtained using method ii) 
above. At any fixed PN order, the normalized region of 
validity w/wisco remains roughly constant as a function 
of a. With a few exceptions, the most conservative es- 
timate vi (lower edge of the error bars in the plots) is 



^ Note that in the Erratum of Paper I we impose a slightly different 
condition; S{vi) = So/2, &{v) = 25o. 



typically in the range v/visco G [0.3 ,0.6]. This is consis- 
tent with the left panel of Fig. O where we see that all 
PN approximations (including high-order ones) peel off 
from the numerical flux in this range. 

These figures allow us to arrive at an interesting con- 
clusion. When we recall that fisco increases with a (cf. 
Fig. [1]), the figures suggest that spin-dependent correc- 
tions in the PN expansion of Eq. ([IJ are effective at push- 
ing the validity of the PN expansion to higher values of 
v/c. However, there is an intrinsic limit to what is achiev- 
able, which is determined instead by w/fisco; and roughly 
independent of a. In the range a € [0.3,0.9], Wisco in- 
creases from ~ 0.444 to ~ 0.609. Therefore the region of 
validity for the orbital velocity is approximately in the 
range v/c € [0.44 x 0.3 , 0.61 x 0.6] ~ [0.13 , 0.37]. 

Let us now focus on the counterrotating case, i.e. on 
the dashed red error bars in Fig. [4l which, again, were 
determined using method ii) above. The only exception 
is the case TV = 3 (corresponding to the right panel of 
Fig. [3]) , that we will discuss separately below. As in the 
corotating case, the region of validity shrinks mildly or 
remains roughly constant as jaj increases. For N = 6 the 
region of validity shrinks faster with increasing spin. 

The edge of the region of validity can also be presented 
in terms of the Boyer-Lindquist radius of the particle's 
circular orbit. The corresponding plots for the corotating 
and counterrotating cases are presented, for complete- 
ness, in Fig.m The ISCO radius rjsco is a monotonically 
decreasing function of the spin (or of i;isco)j so, quite nat- 
urally, the trend as a function of a is the opposite of what 
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FIG. 5. Edge of the region of validity expressed in terms of the Boyer-Lindquist radius for different PN orders, in the corotating 
(black, straight line) and counterrotating (red, dashed line) cases. The blue, dashed lines for A'^ = 3 refer to the counterrotating 
case, and they were obtained by an alternative method (see the discussion around Fig. [7] below). 



we observed for velocities in Fig. 01 Our results consis- 
tently suggest that the region of validity of the PN ap- 
proximation cannot be extended all the way down to the 
ISCO, contrary to a rather common assumption in GW 
data analysis. Instead, one should use care when using, 
for example, the 2PN approximation for r/rjsco < 4, as 
in that regime higher-order PN terms cannot be neglected 
(this is particularly true for rapidly rotating BHs in pro- 
grade orbits). Our results suggest that a safer choice 
would be to truncate all analyses at r/rjsco = 6, which 
ranges between r/M G [6 , 54] depending on the BH spin, 
unless one is dealing with approximants more accurate 
than Taylor expansions. 

Finally, one can also investigate how the edge of the re- 
gion of validity behaves with PN order. This is depicted 
in Fig.inifor a set of fixed values of a (shown with different 
colors, as described in the caption). The vertical dashed 
lines separate the different- iV orders. If we concentrate 
on the nonspinning case (black), ignore the pathological 
case N = 3 (discussed below) and consider the conser- 
vative, lower end of the error bar, we see that there is 
a maximum at iV = 6. For larger values of iV, v/visco 
would consistently decrease, as found in Paper I. In the 
spinning case, however, this trend is not as clear, as at 
N = 6 the edge of the region of validity is rather sen- 
sitive to the spin value. Without higher-order terms in 
the PN expansion, which would provide larger- points 
in this figure, one cannot conclude whether = 6 is the 
optimal order of expansion in the spinning case. 

Before moving on to the next Section, let us discuss 
the = 3 case for counterrotating orbits in more detail. 
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FIG. 6. Edge of the region of validity as a function of A'^ 
for fixed values of a. The thick, black error bar corresponds 
to a = 0, followed by a = 0.3 , 0.6 , 0.9 to the right and a = 
-0.1 , -0.2 , -0.6 , -0.9 , -0.99 to the left. Note that A^ = 3 
is a special case. For reasons discussed in the text, we do not 
give error bars for A'^ = 3 and a < —0.1. 



This is a special case, as noted by the discontinuity in 
counterrotating orbits shown in Figs. 01 El The patholo- 
gies explained below are the reason why, in Fig.jSl we only 
plotted the counterrotating edge of the region of validity 
when the Kerr spin parameter \a\ < 0.1. Notice also that, 
when N = 3, the error regions in Fig. |S1 are significantly 
larger than for any other N value. We will discuss the 
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reason for this below, but the impatient reader can skip 
to the next section without loss of continuity. 

Figure [7] clarifies the origin of the problem. When we 
use method ii) , the top margin of the edge of the region 
of validity is estimated as the (smallest) value of w/wisco 
for which 6'^^''{v) ~ (5o- This condition corresponds to 
the leftmost intersection of the horizontal dashed red line 
with the solid black line in the plot. Similarly, we de- 
termine the most conservative estimate of the edge of 
the region of validity by considering the smallest vi such 
that = 5q/2. This corresponds to the leftmost 

intersection of the horizontal, dot-dashed green line with 
the solid black line. For corotating orbits, as it hap- 
pens, these intersections always exist. In fact, the local 
maximum in 5^'^\v) (which is located at w/wisco ~ 0.8 
for a = 0) moves to the right and becomes significantly 
larger as a — >■ 1. For counterrotating orbits the trend is 
the opposite: the local maximum moves to the left and 
decreases in magnitude. For a critical value of the spin 
a ~ —0.1, the red dashed line and the solid black line 
do not intersect anymore. This is why in Figs. U and [5] 
we only plot the red-dashed (counterrotating) edge of the 
region of validity when \a\ < 0.1. Of course, we can insist 
to identify v and vi as the smallest values of v such that 
6^^\v) = 6o, 5^^^{vi) = (5o/2. This procedure leads to the 
red, dashed error bars in the iV = 3 panel of Figs. 2] and 
[SJ Note that these error bars are unnaturally small for 
\a\ > 0.4. 
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FIG. 7. The solid black line shows S^'^\v) for iV = 3. The 
horizontal dashed red (dot-dashed green) lines represent So, 
computed from Eq. (|10|) . and So/2. See text for discussion. 

Another possible solution would be to switch to 
method i) when method ii) fails. Now the upper mar- 
gin of the edge of the region validity would be given by 
the first zero of S^^'^{v), and the lower margin would be 
estimated by the condition given in Eq. This results 
in the blue, dotted error bars shown in the central top 
panel of Figs. |4] and [5j These error bars are significantly 
more optimistic than the ones we presented in the rest of 
the paper, but (in our opinion) their significance is not 
as clear and well-justified as the rest of our results. 



The problem discussed in this section concerns coun- 
terrotating orbits and TV = 3. This is an exceptional 
case, and it does not affect the conclusions drawn earlier 
in the paper. However we should remark, for complete- 
ness, that similar pathologies occur for corotating orbits 
with = 6 when < a < 0.2. and they may also occur 
at higher FN orders. 



IV. RELEVANCE OF MULTIPOLAR 
COMPONENTS AS A FUNCTION OF SPIN 

Until now, we compared the FN approximation to nu- 
merical results that were considered to be virtually "ex- 
act". This was justified because the Teukolsky code 
computes as many multipoles in the angular decompo- 
sition of the radiation as needed to achieve an accuracy 
of O(10~^°) at any given orbital velocity. While this is 
manageable in frequency-domain calculations, sometimes 
accurate calculations of a large number of multipoles are 
not possible in extreme mass-ratio time-domain codes, or 
in numerical relativity simulations of comparable-mass 
binaries: cf. [13; [13 for an analysis of multipolar decom- 
positions of the radiation from comparable-mass binaries 
and p^ . [20| for more recent numerical work to overcome 
these difficulties. As advocated in several papers [H, 0- 
H, E HH, HI, EMRIs provide a simple playground to 
study the number of multipolar components required to 
reach a given accuracy in the FN approximation (or in 
one of its resummed variants). 

Figure [8] shows a comparison of the convergence of the 
multipolar decomposition versus the convergence of the 
FN expansion of the energy fiux. This plot generalizes 
Fig. 7 of Paper I in two ways: (i) it uses more accurate 
numerical data, and (ii) it considers the effect of the cen- 
tral BH spin on the number of multipolar components 
required to achieve a given accuracy. 

We fix three values of the orbital velocity (v — 0.01, 
i; = 0.1 and v = 0.2) and we plot F^^^' -F^^', where F^^'> 
is the A^th approximant of the FN energy flux and F*^^^ 
is the numerical energy flux truncated at the ^th angular 
multipole. Some features are immediately visible from 
this plot: 

(i) Even at low orbital velocities {v = 0.01), it is nec- 
essary to include multipolar components up to and 
including ^ = 4 to achieve an accuracy better than 
10"'' in the flux; on the other hand, including up to 
^ = 5 we obtain results that arc as accurate as those 
that would be obtained including more multipoles. 

(ii) For an orbital velocity t; = 0.1 (f = 0.2) the 
best-known FN flux and numerical calculations al- 
ways disagree at levels of ~ 10~® (~ 10~^) or 
larger. This is obviously due to the slower, non- 
monotonic convergence of the FN approximation in 
this regime. Some nontrivial features of the FN ap- 
proximation are again well-visible here: for exam- 
ple, as pointed out repeatedly in this paper, when 
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Relevance of multipolar components up to £ = 2, ^ = 3, ^ = 4, ^ = 5, and summing as many ^'s as necessary for the 
accuracy of the Teukolsky code to be 0(10"^°) at any given velocity. 



a = -0.99 and v = 0.1 the 3PN {N = 6) expan- 
sion performs much better than higher-order ex- 
pansions. This may vifell be accidental, and in fact 
it does not hold when v = 0.2, as then iV = 4 is 
(most likely accidentally) better. 

(iii) As a rule of thumb, the inclusion of multipoles up 
to and including £ = 5 seems necessary to achieve 
the level of accuracy of the best-known (A^ ~ 8) 
PN expansion of the flux in the Kerr case. This 
conclusion is independent of the spin parameter a. 
In fact, a has hardly any effect on the number of 
multipolar components that must be included in 
the flux to achieve a desired accuracy. 

V. CONCLUSIONS 

We extended the method proposed in Paper I to de- 
termine the formal region of validity of the PN approx- 
imation for quasicircular EMRIs of compact objects in 
the equatorial plane of a Kerr BH. The boundary of the 



formal region of validity is defined as the orbital veloc- 
ity where the "true" error in the approximation (relative 
to high-accuracy numerical calculations) becomes com- 
parable to the series truncation error (due to neglecting 
higher-order terms in the series). 

For quasicircular, equatorial Kerr EMRIs, the PN ex- 
pansion is known up to 4PN, and our estimate of the 
region of validity can only be pushed up to 3PN. Our 
main results are shown in Fig. 2] (in terms of orbital ve- 
locity) and in Fig. [S] (in terms of orbital radius). At fixed 
but arbitrary spin parameter a, the 3PN approximation 
has no obvious advantage when compared with other PN 
orders. At fixed PN order N, Fig.|3]shows an interesting 
trend: when normalized by the ISCO velocity Wiscoi to a 
very good approximation the region of validity does not 
depend on a. 

We should emphasize that our results say nothing 
about the absolute accuracy of the PN approximation: 
they only suggest relational statements between the A^th 
and the (A^ -|- l)th-order approximations. For velocities 
within the region of validity of the asymptotic series, all 
we can say is that the A^th order approximation has errors 
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that are of expected relative size. For larger velocities, 
the (7V + l)th- and higher-order terms become important, 
and should not be formally neglected. If we can tolerate 
errors larger than those estimated by the (iV + l)th order 
term (at the risk that higher-order approximations may 
be less accurate than lower-order ones) we can surely use 
the PN expansion beyond the realm of its formal region 
of validity. This, however, would force us to lose ana- 
lytic control of the magnitude of the error, as given by 
the next order term. The meaning of this caveat is well 
illustrated by the counterrotating case with a = —0.99: 
it is clear from Fig. [5] that the 3.5PN and 4PN approxi- 
mations do not represent an improvement over the 3PN 
approximation (and in fact perform quite badly) beyond 
the realm of the region of validity. 

Future work could concentrate on studying whether 
resummation techniques, such as Pade [l3 - [l4| or Cheby- 
shev [l^ , enlarge the formal region of validity, using the 
methods developed here. This would allow us to identify 
optimal resummation methods, which, in turn, affects re- 
summed waveform models for EMRIs, like the effective- 
one-body approach [^-Q. Moreover, one could attempt 
to establish whether there is a correlation between the 



accuracy of the energy flux, as measured by the asymp- 
totic methods developed here, and the accuracy of the 
waveform model, as required by GW detectors. 
Note added in proof. After submission of this pa- 
per, we learned that Fujita has extended PN calcula- 
tions of the energy flux up to 14PN for nonrotating 
(Schwarzschild) BHs [2|. 
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